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Structural analysis of branched DNA molecules (BDM) is important as model systems for DNA junctions 
and also as building units for DNA assembly. Although there have been efforts to study the structures of 
BDM, label-free solution structures have not been well determined yet. Here, we used a combination of 
synchrotron-based experimental tools and computational simulation to study the global structures of 
label-free BDM in solution. Overall structures of 3-arm and 4-arm BDM were revealed as an asymmetric 
T(or Y)-shape and a distorted X-shape, respectively. The internal structures of the DNA double helix were 
shown to have a canonical B-form for both the BDM. We also reconstructed the thermal denaturation 
process of BDM by determining the transient global structures over a wide range of temperatures. The 
proposed high-resolution structures of BDM are expected to provide fundamental information for studies 
of the biological function of junction DNAs and DNA assembly. 

Branched DNA molecules (BDM) are widely employed models for studying DNA junctions, which are 
central intermediates in important biological processes, such as genome duplication, DNA damage repair, 
and homologous recombination of DNA''^. In addition to these biological applications, BDM have also 
been used more recently, with the development of DNA technology, as important building blocks in the con- 
struction of complicated DNA nanostructures via molecular self-assembly^ Because the increased importance 
of BDM are due to their structural significance, many studies have aimed to determine the structures of BDM 
using techniques such as atomic force microscopy'^, fluorescence resonance energy transfer (FRET)'""'^, and X- 
ray crystallography" '"'. Despite these efforts, however, determining high-resolution, three-dimensional (3D) 
structures in aqueous solution, where DNA molecules are dispersed, remains a challenge due to technical 
limitations'^". Sample immobilization, trapping, chemical labeling, or crystallization, which are required for 
using these structural tools, prohibit the discovery of realistic structures of branched DNA. Structures determined 
by these methods allow only limited information because of a lack of freedom in configurations of DNA 
molecules, especially in the aqueous phase. 

Solution X-ray scattering is a powerful method to determine solution-phase, label-free structures of macro- 
molecules. The structure determination based on the assumption of single molecular envelop has become 
increasingly convenient owing to the compilation of shape determination tools"* '". In addition, DNA assembly 
is generally performed in solution phase. However, due to the intrinsically reduced information, which is one- 
dimensional and orientationally averaged scattering intensity profile, difficulties still remain in probing systems 
with conformational flexibility and multiplicity^". To address this point, ensemble-based reconstruction methods 
have been suggested that offer insight into the conformational transition^'"^'*. 

In this report, we show that the combination of synchrotron smaU-angle X-ray scattering (SAXS), wide-angle 
X-ray diffraction (WAXD), and atomistic simulations can be powerful in structural studies of label-free, BDM 
dispersed in aqueous solution. By using synthetic 3-arm and 4-arm BDM as model systems for three-way 
junctions (3 WJ) and four way junctions (4 WJ), respectively, we have reconstructed the global structures of 
junction DNAs in the aqueous phase. In addition to determining high-resolution structures at room temperature, 
we have also determined the transient global structures of BDM over a wide range of temperatures to study the 
denaturation process of 3 WJ and 4 WJ DNA. 
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Results 

In order to synthesize monodisperse BDM, we employed the one-pot 
synthesis method'"'^: equal moles of all 3 and 4 oligonucleotides were 
mixed together to form the 3-arm and 4-arm BDM (hereafter, 
referred to as 3 WJ and 4 WJ DNA, respectively). The formation 
of BDM was evaluated by gel electrophoresis (Supplementary Fig. 
SI), in which the mobility of DNA molecules depends on their size, 
shape, and extent of base pairing. The estimated yield of branched 
DNA is close to 100%. The sequences of the oligonucleotides are 
presented in Supplementary Table SI. 

Determination of 3 WJ and 4 WJ DNA structures in aqueous buf- 
fer at room temperature by solution SAXS. Synchrotron SAXS 
measurements of 3 WJ and 4 WJ DNA in PBS buffer were carried 
out at 25^C (Fig. 1). Guinier analysis^'', performed on the SAXS data 
(Supplementary Fig. S2), revealed linearity in a small q region which 
suggests that the DNA molecules exist in a monodisperse conforma- 
tional state at each temperature. In Figure la, we display SAXS 
intensity curves, I(q), measured for 3 WJ and 4 WJ DNA at 25°C. 
To provide an intuition for overall molecular shape in real space, the 
pair-distance distributions, p(r), are shown in Figure lb, which were 
calculated by the indirect Fourier transform of the entire scattering 
profiles using the program GNOM^' To further facilitate direct 
structural comparison, we created atomistic structures of 3 WJ and 
4 WJ DNA by molecular dynamic (MD) simulations. The aforemen- 
tioned quantities, I(q) and p(r), can be computed from the atomic 
coordinates as well, through the use of the program CRYSOL^" and 
GNOM^', respectively. The calculated SAXS curves and pair- 
distance distributions for representative configurations of 3 WJ 
and 4 WJ DNA obtained from MD simulations are generally in 
agreement with the experimental results at 25°C, as shown in 
Figure la and lb. For a quantitative comparison for p(r), we ob- 
tain the radius of gyration, Rg,p(r)) and the maximum dimension, 
Dmaxj as the distance where the p(r) function approaches zero. The 
results for Rg,p(r) and D^ax obtained from experiment at room 
temperature agree reasonably with those from the simulation 
results, as listed in Supplementary Table S2. We observe that both 
the experimental and simulation results of p(r) exhibit three main 
peaks characterizing geometry of BDM in Figure lb. The first peak at 
r ~ 2 nm reflects average distances between backbones, correspond- 
ing to the average diameter of double helix. The second one around r 
~ 3.5 nm, although not distinct, attributes to atom pairs constituting 
the pitch of helix turn. The third peak located at r ~ 5 nm 
corresponds to average distance between each helix arms, which 
determines the overall size of molecular envelope. 

To corroborate the p(r) function analysis and to obtain a visual 3D 
representation of the structure of 3 WJ and 4 WJ DNA in different 
conformational states, we reconstructed model-independent struc- 
tures that occupy a search volume defined by the maximum dimen- 
sion, Dnijix, calculated from the p(r) function using the ab initio shape 
determination program DAMMIF^', an upgraded version of 
DAMMIN™. Multiple runs of DAMMIF were used to generate struc- 
tural models, with and without imposing symmetry constraints (No 
symmetry, P2 symmetry, P3 symmetry, and P4 symmetry). The 
structural models for 3 WJ DNA in PBS buffer at 25°C exhibited 
distorted Y-shaped forms composed of 3 arms and the junction in the 
innermost center of the reconstructed models, except in that with 
imposed P3 symmetry restriction (Fig. Ic). The structural models for 
4 WJ DNA in PBS buffer at 25°C adopted distorted X shaped forms 
consisting of 4 arms and the junction near the center of the recon- 
structed models, which slightly deviated from an ideal X-shaped 
conformation with imposing P4 symmetry constraint. Some struc- 
tural differences in the overall shape were observed depending on the 
type of symmetry restriction imposed (Fig. Id). The reconstructed 
conformations resemble the atomic structures obtained from MD 
simulation for both 3 WJ and 4 WJ DNA in Figures le and If It is 



fascinating that 3D structures of DNA molecules in solution were 
successfully determined from X-ray scattering experiments without 
any constraint such as labeling or trapping. 

A detailed description of equilibrium structures obtained from 
MD simulations is provided in Supplementary Information (Figs. 
S3 and S4, and Table S3). In particular, we carefully compared the 
results obtained from 2 force field parameters that have been widely 
used. We also discuss the discrepancies between the experimental 
results and simulation results in the p(r) patterns in detail, by means 
of the relationship between the general features in the p(r) function 
and the overall geometry of the DNA molecules (Figs. S5 and S6). 

General features of branched DNA by solution WAXD. With the 
aim of better understanding the helical structures within DNA 
molecules, synchrotron WAXD measurements of 3 WJ and 4 WJ 
DNA in PBS buffer were performed at 25°C (Fig. 2). To gain fur- 
ther insight into the internal helical structures of DNA molecules, we 
evaluated ensemble-averaged WAXD patterns of the 2,500 model 
configurations calculated from MD simulations using the program 
CRYSOL. The calculated WAXD patterns can be compared with the 
experimental WAXD patterns at 25"C by means of the positions and 
intensities of diffraction peaks, labeled P1-P6 in Figure 2. Scattering 
from structural repeats in helical DNA molecules produces a series of 
diffraction peaks in the WAXD pattern that vary in position and 
intensity, depending on the molecular structure" ''^. Figure 2 shows 
ensemble-averaged WAXD patterns of full (structure 3 in Fig. 2) and 
partial (structures 4 and 5 in Fig. 2) structures of the 2,500 model 
configurations of a 3 WJ DNA molecule calculated from MD 
simulations. Scattering caused by the stacked nucleotide bases can 
be estimated from the partial structure 4 in Figure 2, which yields a 
single diffraction peak, P6, shown by the orange line. Peaks in the 
small angle region, P1-P4, originate from scattering due to the sugar- 
phosphate backbone, as shown by the green line in the WAXD 
pattern calculated from the partial structure 5 in Figure 2. How- 
ever, we do not observe a diffraction peak, P5, that arise predomi- 
nately from cross-terms between the backbone and stacked base 
pairs which are enforce in a B'-form DNA conformation'''. As 
shown in Figure 2, the experimental WAXD patterns for both 
3 WJ and 4 WJ DNA molecules in PBS buffer at 25°C exhibited 
P1-P6 diffraction peaks, with the exception of P5, that are charac- 
teristic of a typical B-form DNA conformation. 

Conformational change of branched DNA in thermal denatura- 
tion. To fuUy understand the conformational behavior of DNA, it is 
necessary to define the structure of both the native state and the 
various denatured states with the change of environments. Recent 
studies have revealed the biological significance of denatured states, 
functional flexibility, and intrinsic disorder that occurs in functional 
regions and interfaces under physiologically active conditions'^''. 
However, the detailed structural features of BDM and the varia- 
tions in structure that occur in response to temperature change 
have not been well studied. 

To gain a visual 3D representation of the overall structures and the 
structural differences of 3 WJ and 4 WJ DNA molecules in solution 
at various temperatures, we attempted to reconstruct structural mod- 
els using the ab initio shape determination program DAMMIF. The 
experimental SAXS curves and their p(r) functions at temperature 
range between 25°C and 72°C are shown in Figure 3. As shown in 
Figures 3a and 3b, the scattering profiles changed in the wide q region 
(q > 1.5 nm"') between 63°C and 69^C. The scattering patterns at 
72^C were clearly changed in the whole q region, due to conforma- 
tional transition from the branched DNA duplex structure to a fully 
random coil conformation of single-stranded DNA. These results 
were consistent with the disappearance of the 3rd peak, which is 
characteristic of a branched DNA form, in the p(r) function (Supple- 
mentary Table S2). The structural parameters obtained from the 
SAXS data of 3 WJ and 4 WJ DNA molecules in solution at various 
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Distorted Y Shape ^0A_ Distorted X Shape 

Figure 1 | SAXS data and structural models, (a) Scattering profiles of 3 WJ and 4 WJ in PBS buffer at 25°C. Insets show 2D SAXS patterns. The symbols 
are the experimental data, and the solid lines are SAXS profiles obtained from the reconstructed structural models using the ab initio shape determination 
program DAMMIF. The dashed lines without symbols are the theoretical SAXS curves of structural models calculated from MD simulations. For clarity, 
each plot is shifted along the log I axis, (b) Pair distance distribution functions for 3 WJ and 4 WJ, based on an analysis of the experimental SAXS data 
using the program GNOM. The dotted lines without symbols are the p(r) functions for structural models obtained from MD simulations labeled as 
dashed lines in a. (c-d) Reconstructed structural models of 3 WJ and 4 WJ in PBS buffer with and without any symmetry restrictions from the ab initio 
shape determination program DAMMIF DAMMIF (Chi = 1.599 (PI), 1.603 (P2), 1.611 (P3) for 3 WJ and Chi = 1.545 (PI), 1.549 (P2), 1.544 (P4) for 
4 WJ). (e-f) Structural models calculated from MD simulations for 3 WJ and 4 WJ using simulation employing the parmbscO force field. Surface 
rendering in the structural models was achieved using Discovery Studio 1.6 (Accelrys Inc.). 

temperatures are listed in Supplementary Table S2. On the whole, the a substantial conformational change at 72°C during the thermal 
obtained values of both DNA molecules at 72°C were markedly denaturation process. Such observations which occurred in the ther- 
decreased in comparison with those measured at 25-69 C. From mal process could be explained in detail by the 3D reconstructed 
these results, we observed that both 3 WJ and 4 WJ DNA underwent models. Figure Be, f show the structural models reconstructed for 



SCIENTIFIC REPORTS | 3 : 3226 | DOI: 10.1038/srep03226 



3 



d (A ') = 2ji/q 
12.5 6.3 42 3.1 2.5 



P1 P2 P3 P4 P6 




"T 1 1 ' 1 ' 1 ' r— 

5 10 IS 20 25 

q (nm ') 




3 4 5 

Figure 2 | WAXD profiles from experimental studies and MD 
simulation. WAXD fingerprint patterns from 3 WJ (blue) and 4 WJ (red) 
in PBS buffer at 25°C. Ensemble-averaged WAXD patterns calculated from 
MD simulations for the complete structure of 3 WJ (black), the base pair 
atoms only (orange), and the sugar-phosphate backbone atoms only 
(green). 

3 WJ and 4 WJ DNA molecules in PBS buffer at various tempera- 
tures using D AMMIF. The reconstructed models were obtained with 
and without imposing any restrictions on the symmetry of the DNA 
molecule. As shown in Figure 3, structural differences in the overall 
shape were apparent depending on the symmetry restriction, but all 
models are very meaningful for structure interpretation. However, in 
order to strike a balance between symmetry underestimate from PI 
(no symmetry) constraint and symmetry overestimate from P3 (for 
3 WJ) and P4 (for 4 WJ) constraints, we focused primarily on the 
structural models reconstructed with imposing P2 symmetry restric- 
tion that consistently represent the most reasonable and adequate 
feature for a systematic understanding of the conformational change 
in response to change in temperature. As already shown, the struc- 
tural models reconstructed for 3 WJ DNA in solution at 25°C 
include 3 distinct arms of the same length. However, in the 3 WJ 
DNA structures modeled at 40-63°C, under early stages of the dena- 
turation process, 2 of the 3 arms were somewhat shortened compared 
with those of the 3 WJ DNA at 25°C. These results suggest that the 
3 WJ DNA underwent a slight conformational change during ther- 
mal denaturation, which is rearrangement of the denatured strands 
that was caused by breakdown of hydrogen bonding between bases 
within individual arm domains. In addition, it can be predicted that 
at least 1 arm domain is likely to retain a certain degree of structural 
stability in the overall shape under the present conditions. The struc- 
tural models reconstructed for 3 WJ DNA in solution at both 66 and 
69°C revealed bent forms that were different from those of 3 WJ 
DNA in solution at lower temperature; one of the 3 arms completely 
disappeared due to destruction of the DNA duplex structure. These 
results seem to be correlated to the change of ratio of the 1st peak 
intensity values to the 3rd peak ones in the p(r) function. We 
observed that the contribution of the 3rd peak was gradually 



decreased as the thermal denaturation was proceeding. The 3 WJ 
DNA molecule in solution at 72°C, during the final stages of 
denaturation, adopted a fully random coil conformation of single- 
stranded DNA, which is confirmed by analysis of p(r) for single- 
stranded DNA obtained from MD simulation (Supplementary 
Figs. S7 and S8, and Table S4). The 4 WJ DNA in solution at 25°C 
was found to exist in a distorted X- shaped conformation containing 
four arms of the same size, and was very similar to the simulated 
conformation. The structural models reconstructed for 4 WJ DNA 
in solution at 40 to 60°C, in the early stages of denaturation, have 
similar four-armed conformations overall to those of 4 WJ DNA at 
25°C, but more distorted X-shaped forms. Interestingly, from the 
4 WJ DNA structures modeled at around 63°C, we found that the 
4 WJ DNA experienced a considerable structural change, where the 
molecule changed from a 4-armed to a 3-armed conformation. The 
structural models reconstructed for 4 WJ DNA maintained their 
structure based on the 3-armed conformation up to 69°C. These 
results suggest the possibility that one or two of the 4 arms are 
denatured in advance and then they are aggregated to form a partially 
molten globule. Interestingly, structural models (at 63 to 69'C) 
reconstructed by imposing a P4 symmetry constraint have some 
similarities to those (3-armed conformation) reconstructed based 
on P2 symmetry restriction, while retaining P4 symmetry character- 
istics in the right model of each column. The 4 WJ DNA molecule in 
solution at 72°C, during final stages of thermal denaturation, 
adopted a fuUy random coil conformation of single stranded DNA, 
similar to the 3 WJ DNA in solution at 72°C. From these reconstruc- 
tions, we were able to identify an intermediate DNA conformation 
that was dominant at each temperature in the thermal denaturation 
process. 

In order to understand these thermal denaturation process, we 
additionally performed ultra violet (UV) absorbance and circular 
dichroism (CD) spectra analysis. The thermal denaturation of double 
stranded DNA is progressive and this can be seen in a sigmoidal 
increasing of UV absorption. By increasing temperature, the trans- 
ition from the double stranded to denatured structure was observed 
over the range of approximately 60 to 70°C. The T^j values deter- 
mined from UV absorption spectra are 64 and 63 "C for 3 WJ and 
4 WJ, respectively (Fig. S9). After 72°C, both EDMs are fully dena- 
tured to three and four single stranded oligonucleotides. These ther- 
mal behaviors were also shown from the CD spectroscopic analysis. 
(Fig. SIO). These UV and CD results are consistently well matched 
with previously described SAXS characterization. 

MD simulation results at high temperatures. We have performed 
MD simulation of the atomistic model of BDM at elevated tempera- 
ture in order to obtain a microscopic description of structural 
transitions that occur during thermal denaturation. Our model 
system for BDM solvated in PBS buffer remains stable rigidly for 
longer than 500 ns at 400 K (Figs. SI 1 and S12). In order to facilitate 
denaturation, we increased the temperature of the system equili- 
brated at 400 K to 450 K instantly, which eventually causes melt- 
ing within 100 ns. We note that the temperature applied in our 
model system does not scale consistently with the physical tempera- 
ture. The use of unphysical temperature was intended to apply an 
enhanced noise level to the system. We discuss the simulation result 
at high temperatures in Supplementary Information in detail. 

In Figure 4, the melting process is depicted with isosurface plots 
for the probability density of finding atoms in 3D space, computed 
from three subsequent trajectories of 30 ns for 3 WJ and 4 WJ sys- 
tems. In the experiment, averaging over numerous molecules and 
time is reflected in the diffraction pattern, such that mobile portions 
produce relatively low intensity. Likewise, portions of molecules 
that are relatively unstable should yield a lower probability density 
of atom existence in our analysis. In this manner, the isosurface 
plot is considered a reasonable counterpart of MD simulation to 
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Figure 3 | Structural models in response to change in temperature based 
on small-angle X-ray scattering. Scattering profiles of 3 WJ (a) and 4 WJ 
(b) in PBS buffer as a function of temperature. The symbols are the 
experimental data, and the solid lines are SAXS profiles obtained from 
reconstructed structural models by using the ab initio shape determination 
program DAMMIF. For clarity, each plot is shifted along the log I axis. Pair 
distance distribution functions for 3 WJ (c) and 4 WJ (d) in solutions of 
various temperatures, based on an analysis of the experimental SAXS data 
using the program GNOM. Reconstructed structural models of 3 WJ (e) 
and 4 WJ (f) in PBS buffer, with and without any symmetry restrictions, as 
a function of temperature using the ab initio shape determination program 
DAMMIF (Chi = 1.477 ~ 1.635 for 3 WJ and Chi = 1.447 ~ 1.789 for 
4 WJ). Surface rendering in the structural models was achieved using 
Discovery Studio 1.6 (Accelrys Inc.). 

understand reconstructed structures from the experiment. In 
Figure 4, we additionally show the probability density of finding 
Na* ions. Specifically, the isosurfaces correspond to the position 
where the chance of finding DNA and Na^ ions is 30% and 1%, 
respectively. Pair-distance distributions, p(r), computed from 
instantaneous structures are also presented indicating partially dis- 
ordered structures of 3 WJ and 4 WJ DNA in Figure 4. While the 
molecules undergo thermal denaturation, changed are the positions 
of the peaks in p(r) and size of the molecules. In describing the 
melting process, we focused on 3 features that are common to both 

3 WJ and 4 WJ. After increasing the temperature, enhanced thermal 
noise causes large fluctuations of both local and global structures, 
including inter-helix angles. As a result, a distorted structure of the 
backbone facilitates base opening for 0-30 ns. Interestingly, base 
pairs of 2 helices start to flip out, while those of the other helix remain 
intact. Next, for 30-60 ns, the 2 helices with broken hydrogen bonds 
are shortened to form a partially molten globule. Finally, after 60 ns, 
2 helices are aggregated, leading to the complete collapse of 1 helix. 
This aggregation is not considered the result of any targeted inter- 
action, such as the formation of hydrogen bonds between 2 helices, 
but the cooperative effect of surrounding counterions and water 
molecules. The attractive electrostatic interaction between negatively 
charged phosphate and Na* ions compete with the hydration force 
exerted on Na* ions''*. With the increased thermal noise, the hydra- 
tion effect of water molecules for both the DNA coil and ions is 
reduced, resulting in a higher probability that Na* ions wiU be 
located near DNA, as shown in Figure 4. This observation provides 
a microscopic picture of reconstructed model structures presented in 
Figure 3, especially for shortened or missing helix at the intermediate 
stages of denaturation. 

Discussion 

In summary, we show the entire process of synthesizing and struc- 
turally analyzing solvated 3 WJ and 4 WJ by combining SAXS, 
WAXD experiments, and atomistic simulations over a wide range 
of the temperature. The global structures of BDM are determined 
successfully by interpreting the SAXS intensity profiles through the 
use of shape determination software. Revealed structures are an 
asymmetric T- (or Y-) shape and a distorted X-shape for 3 WJ and 

4 WJ DNA, respectively, at room temperature. We also determined 
from the WAXD experiments that the double helix structures in both 
branched DNAs belong to a canonical B-form. Moreover, our results 
provide an intermediate DNA conformation, which was dominant at 
each temperature during thermal denaturation, by analyzing the 
results from in-situ temperature SAXS in solution and MD simula- 
tions. In the latter, we offer a microscopic insight for thermal dena- 
turation of BDM, taking the cooperative dynamics of counterions 
and water molecules into consideration. The methods applied in this 
article could be used for studying structures of other nucleic acids, 
such as DNA hairpins, enzymes, bent sequences, and RNA, as well as 
proteins and protein complexes. In general, solution SAXS can 
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Figure 4 | Global structures of 3 WJ and 4 WJ derived from MD simulations. Three stages of melting processes in 3 WJ and 4 WJ systems, obtained 
from nonequilibrium MD trajectories, where the temperature is increased from 400 K to 450 K instantly at t — 0. (a) Probability isosurfaces where 
the chance of finding atoms for 3 WJ DNA (blue) and Na"^ ions (magenta), computed with 3 subsequent trajectory segments of 30 ns long, are 30% and 
1%, respectively, (b) Probability isosurface for 4 WJ DNA using the same criteria as in a. In each, the right panels are pair-distance distributions, p(r), for 
individual configuration of 3 WJ and 4 WJ at each stage. 



provide solution structural information at resolutions often suf- 
ficient for functional insights into how these bio molecules work in 
the context of conformational and assembly state under physio- 
logical conditions. The results reported here will guide the develop- 
ment of new DNA nanostructural and self-assembly techniques 
using BDM as building units. 

Methods 

Sample preparation. The sequences of oligonucleotides for BDM were modified 
from Lug's X- and Y-DNA^^-^^. Oligonucleotides were resuspended in TE buffer. The 
3 WJ and 4 WJ DNA were constructed by mixing three or four 36-bp 
oligonucleotides, respectively, in PBS buffer. 3 WJ and 4 WJ DNA were synthesized 
and purified following the methods described previously^-^^'^^. 

X-ray data acquisition and processing. SAXS measurements were carried out using 
the 4 C SAXS II beamline of the Pohang Light Source II (PLS II) with 3 GeV power at 
the Pohang University of Science and Technology (POSTECH), Korea. A light source 
from an In-vacuum Undulator 20 (IVU 20: 1.4 m length, 20 mm period) of the PLS II 
storage ring was focused with a vertical focusing toroidal mirror coated with rhodium 
and monochromatized with a Si (1 11) double crystal monochromator, yielding an X- 
ray beam wavelength of 0.675 A. The X-ray beam size at the sample stage was 0.2 (V) 
X 0.6 (H) mm^. A two-dimensional (2D) charge -coupled detector (Mar USA, Inc.) 
was employed. A sample- to -detector distance of 2.00 m for SAXS and 0.25 m for 
WAXD were used. The magnitude of scattering vectors, q — {4:n/?.) sin 9, was 
0.2 nm~' < q<25 nm~\ where 20 is the scattering angle and a is the wavelength of 
the X-ray beam source. The scattering angle was calibrated with polyethylene-/?- 
polybutadiene-/7-polystyrene block copolymer and silver behenate standards. We 
used solution sample cells with 10 )j.mthickmica windows, a volume of 50 |j.L, and an 
X-ray beam path length of 0.8 mm. AH scattering measurements were carried out at 
various temperatures by using a JULABO FP50-HL control system with a 
temperature stability of ± O.OL C. The SAXS data were collected in 5 successive 
frames of 0.1 min each to monitor radiation damage. There were no changes in the 
scattering patterns with time, i.e., no radiation damage was detected during the 
scattering measurements. To obtain high-quality scattering data without any 
interferences of DNA molecules (ie., concentration effects) (Supplementary Table 
S2), measurements of DNA solutions were carried out over a small range of 
oligonucleotide concentrations (3.2-6.7 mg/mL). In addition, the molecular mass of 



each sample has been calculated to confirm the consistency between the expected 
molecular masses from the sequence of the nucleotides and those calculated from the 
forward scattering 1(0) (in the absolute scale after normalization against 
concentration and water scattering) from Guinier extrapolation (Supplementary 
Table S5). Each 2D SAXS/WAXD pattern was circularly averaged from the beam 
center and normalized to the transmitted X-ray beam intensity, which was monitored 
with a scintillation counter placed behind the sample. The scattering of the PBS buffer 
solution was used as the experimental background. 

MD simulation and structural analyses. We built atomistic structures of 3 WJ and 
4 WJ DNA, corresponding to the sequences displayed in Supplementary Table Si, 
using the NAB tool included in the AMBERtoolsl2 suite'^. Specifically, the structures 
of double helices for each arm were constructed separately in a generic B-form using 
NAB, and arranged appropriately to connect the inter-helical bonds and form the Y- 
shape and X-shape of 3 WJ and 4 WJ, respectively. We then minimized the molecule 
for 2 ps using the steepest descent algorithm, and raised the temperature to 300 K 
gradually over 100 ps to relax the inter-helical bonds, which were initially very 
awkward, while monitoring the bond potential energy. Although not perfectly 
relaxed, the resulting structures of DNA molecules were found to be stable 
throughout the subsequent procedures of solvation and equilibration. 

Most of the present MD simulations were performed using the GROMACS soft- 
ware package^^ with the all-atom parmbscO force fields^^, which is known to preserve 
the stable double helix structure during long simulations^^. For comparison, we also 
carried out simulations using the parm99 force field*". The DNA molecule was 
neutralized with 0.15 MNa^ and CI" ions, and solvated with more than 28000 TIP3P 
water molecules*' in a rectangular box. The size of the simulation cell was set to ensure 
a solvation shell of 10 A thickness around the DNA molecule. We carried out MD 
simulations at constant temperature (300-450 K) and pressure (1 bar) using the 
Berendsen algorithm*^, with a time constant of 1.0 ps for heat bath coupling and 
pressure relaxation. The compressibility was set to 4 X 10"^ bar"\ corresponding to 
the water system. An integration time step of 1 fs was selected, and the electrostatic 
interactions were calculated using the Particle Mesh Ewald (PME) method*^ with a 
real space cut-off of 9 A. Initially, the DNA molecule was relaxed again through the 
energy minimization with water molecules and ions, and allowed to equilibrate at 
300 K at constant volume for 500 ps. Subsequently, the equilibration was continued 
at constant pressure for at least 100 ns, whUe fluctuations of the root-mean square 
deviation were monitored. The trajectories obtained from production runs thereafter 
were used for analyses. 
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